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Abstract 

This is the third paper of a series revisiting the Faraday effect. The question of the absolute 
[ convergence of the sums over the band indices entering the Verdet constant is considered. In 

^Sj ■ general, sum rules and traces per unit volume play an important role in solid state physics, 

and they give rise to certain convergence problems widely ignored by physicists. We give a 
j^', complete answer in the case of smooth potentials and formulate an open problem related to 

less regular perturbations. 
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1 Introduction 



In this paper, vifhich is the third in a series, we continue our study begun in |Bj where we outlined 
a general and mathematically sound theory of the Faraday effect. The theoretical investigation of 
this effect -which amounts to the study of the transverse conductivity of the optical sample- has 
' a long and distinguished history in the quantum theory of solids. While sending the reader to [B] 

. for a detailed discussion and more references, we only mention here that the first general theory 

of the conductivity tensor for independent electrons subjected to a periodic electric potential 
and to a constant magnetic field in the linear response approximation has been developed by L. 
Roth [TB]. Roth's method is based on an effective Hamiltonian approach for Bloch electrons in 
a weak magnetic field (see [T7] and references therein). But due to the difficulties generated by 
QQ ' the linear growth at infinity of the magnetic vector potential, the theory is highly formal. Basic 

, questions as for example computation of higher order corrections in the magnetic field strength, 

the thermodynamic limit as well as some convergence issues, seem to be almost hopeless to cope 
with. Accordingly, the obtained formulae -although considered as a landmark of the subject- still 
' contain singular terms at the crossings of the Bloch bands and only met a moderate success at 

, the practical level. 

In the first paper of this series we outlined a rigorous theory of the Faraday effect free of all these 
difficulties. The first basic idea is to employ a method going back at least to Sondheimer and Wilson 
[T8j in order to express the traces involved in computing various physical quantities as integrals 
using Green's functions, which are easier to control than the eigenfunctions and eigenvalues. The 
^ ' second and the main tool is the use of gauge covariance via a regularized magnetic perturbation 

5^ ! theory [Il[3lll[5l[7l[i[l0l[ni[Il[l3]. The result is a well defined, gauge invariant, workable 

formula for the transverse conductivity tensor in the thermodynamic limit. This tensor has an 
asymptotic expansion in the strength of the magnetic field, whose first two coefficients have been 
explicitly written down in terms of the zero field Green function, which in turn can be expressed 
in terms of zero field Bloch energies and eigenfunctions (see Section VII in [B]). 

The subsequent papers of our series are devoted to some more delicate mathematical aspects 
of the program left open in [6]. In particular, in the second paper of the series [7] we solved -by 
making use of the full power of the regularized magnetic perturbation theory- the most difficult 
part of the program: the proof of the existence of the thermodynamic limit. 
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The present paper is devoted to another convergence issue left open in [B] . More precisely, when 
expressing the relevant quantities in terms of Bloch eigenfunctions and eigenvalues one arrives, 
by formal computations, at expressions containing multiple sums over the band indices as well as 
integrals over Brillouin and unit cells. The point is that the order of summation is crucial since 
the sums might not be absolutely convergent. The problem is more than an academic one since in 
practical computations finite bands models are used which amounts to truncate the infinite sums 
so that the question whether the sums are absolutely convergent is a crucial one. The main result 
of this paper is that the smoothness of the periodic potential V guarantees that such series are 
always absolutely convergent. The smoothness of V is not superfluous as shown by an example. 
Since the same kind of series appear very oftenly in solid state physics (e.g. when computing 
derivatives of dispersion laws) and are related to sum rules concerning various observables, we 
present our main technical result Theorem 13. II in a general setting . 

The content of the paper is as follows. Section 2 contains the setting, examples of series we are 
interested in as well as the related sum rules. In Section 3 we give our main technical result giving 
the decay (under suitable smoothness condition of V) of the matrix elements of the momentum 
operator in the Bloch basis as the band index increases. Section 4 contains applications to series 
appearing in the theory of the Faraday effect. Finally in Section 5 we present a computation 
involving the (singular) one dimensional 5-" potential" , in order to point out the need of smoothness 
conditions on V. The question to optimally relate the decay properties of the momentum matrix 
elements in the Bloch basis to the degree of smoothness of V is left as an interesting open problem. 



2 Preliminaries and motivation 

The one particle Hilbert space for a d-dimensional spinless particle is Hoo = L'^{M.'^), d > 1. The 
Hamiltonian describing its dynamics is: 

H^^P^ + V, P = -iV, (2.1) 

where we assume that V is e- relatively bounded with respect to — A, and periodic with respect to 
a Bravais lattice C. Denote by fl its elementary cell and by fl* the corresponding Brillouin zone. 

We assume for simplicity that the lattice is Z'', rf > 1, and the unit cell il is the unit cube 
centred at the origin. The Bloch-Floquet unitary is (see e.g. |15|): 

U : L^{R'^) / L'^{n)dk, 
Jn* 

(f//)(x,k):=— 1^ ^ e-''-(^+^)/(x + 7), kef7*, x^n, / e C„°°(R'^). (2.2) 

Then H^o can be written as a fiber integral H^o = J^, /i(k)(ik where the fiber operator 

h{k) = (-zVp + k)^ + V, (2.3) 

is defined in L^(fl) with periodic boundary conditions, and its domain is the periodic Sobolev 
space H'^{flp). Here —i'Vp denotes the momentum operator in L'^{fl) with periodic boundary 
conditions. The operator h(k) has compact resolvent and purely discrete spectrum {Aj(k)}j>i. 
We can choose a set of eigenfunctions {Mj(-, k)}j>i which form an orthonormal basis of L^(f7) and 
obey the equation: 

Mk)^.,(.,k) = A,(k)^,(.,k). (2.4) 

We label Aj(k) in increasing order. We have to remember that due to this choice, Aj(k) and 
Mj(x, k) are not differentiable with respect to k at crossing points. Without loss of generality, we 
may also assume that /i(k) > 1 for all k which implies: 

Aj(k)>l, ken*, j>l. (2.5) 
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One can introduce the Bloch functions: 

M/,(x,k) := -^«,(x,k), xeK'* (2.6) 

which form a basis of generahzed eigenfunctions of ffoo in L^{^'^), i.e. in distributional sense we 
have: 



/ ^M/,(x,k)vI/,.(y,k)dk = 5(x-y). 



Using Bloch functions, one can express the integral kernel of any function of H^o in terms of the 
fiber /i(k). For example, the Green function (i.e. the integral kernel of the resolvent) writes as: 



G„,x.y;,)=/ S2i(:Jg^<ik, (2.7) 

The above formula has to be understood in the distributional sense since the series on the right 
hand side is typically not absolutely convergent. Now let us present the problems we want to 
consider. 



2.1 Perturbation theory and sum rules 

By writing 

/i(k) = /i(ko) + 2(k - ko) • [-iVp + ko] + (k - ko)2 =: /i(ko) + W, 

we see that W is a regular perturbation for /i(ko). Assume that Ai(ko) is non-degenerate. Then 
according to the analytic perturbation theory, Ai(k) remains non-degenerate in a neighborhood 
of ko. Keeping |k — ko| small enough and using the Feshbach formula with the projection Ho = 
|'Ui(-, ko))(ui(-, ko)| we obtain: 

A(k) = A(ko) + (ui(-,ko),VFui(-,ko)) 

- {Wui{; ko), {n^m - A(k)]n^}"' Wui{; ko)). (2.8) 

By iterating the above formula we can identify in an efficient way the full Taylor expansion of Ai 
around ko. Let us first introduce the notation: 



7rjj(a,k) := (Mj(-,k), (-i^a + fca)Mj(-,k))i2(n) = / Ui{x,k)[{-ida + ka)uj]{x,k)dx. (2.9) 

Jn 

If we are interested for example in ^p^, one of the terms building it will be proportional with: 

2. -J-Tx^ 2^ 7r,3i(a,k), (2.10) 

where the order of the sums is crucial. Each series generates an Z^(N)-summable vector for the 
next one. The multiple series is convergent, but not apriori absolutely convergent. This is because 
the inner tt's can grow in absolute value with the energy. But we will prove in this paper that 
such series arc always absolutely convergent if V is smooth. 

Before giving some more precise results, let us make the connection with the notion of sum 
rules and argue why the above matrix eleincnts may have a rapid decay with the energy if one 
index is kept fixed. We work in L^(ri) and let us omit k in order to simplify notation. If A is an 
operator which is relatively bounded to h, then we can write the identity {t> 0): 

i([Ae'*''^e-'*^]u„,n„) =2^|(^u„,w„)|2sin{t(A„-A„)}. (2.11) 

n>l 
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The commutator on the left hand side is well defined as a quadratic form on functions belonging 
to the domain of h. Now if we assume that hAh~^ is bounded for some > 1, then we obtain 
the bound 



which when applied to (j2.11D would insure that the series defines a differentiable function at all t. 
If in addition we assume that we can make sense out of the double commutator [A, [/i, A]], then a 
typical sum rule would be: 



([A [h^A]] 

Urn 1 '^m/ — 2 \{Au,n, ■ (An - Am). (2-12) 

n>l 



If A does not commute well with h, it can happen that even if the function of t in (j2.1ip is 
continuous on R, we cannot be sure that it is also differentiable. We will later on give an example 
of a Schrodinger operator which will make appear the Riemann function on the right hand side 
of (|2.12p . This function is known to be everywhere continuous, but differentiable only at certain 
rational points not including t — Q. 

For many more rigorous aspects and applications of sum rules, see [5]. 

2.2 Traces per unit volume 

Let us consider the Fermi-Dirac distribution function: 

fpDix) := ^^(,_],)^^ , a;elR,/3>0,Ai€]R, (2.13) 
and choose the following path in the complex plane 

T = {x±i5: -1 < X < cx)} IJ {-1 + iy: -5 <y<5} (2.14) 
with a fixed 6 £ (0, Introduce the notation: 

Ro^iz) := (H^ - z)-\ (2.15) 
The quantity we are interested in this time is a trace-per-unit- volume (provided it exists): 

:= 1^ y rfx |y dz } FD{z)PatRoo{z)Pa^Roo{z) . . . Pa,,Roo{z)^ (x, x). (2.16) 

We want to express the above quantity only with the help of Bloch functions and energies. The 
fact that the operator 

/ {-^Qi -^oo (^)-^a2 

Roo{z)...Pc.^Roo{z)} fFD{z)dz (2.17) 

has a jointly continuous integral kernel if V is smooth, can be proved with the same methods as 
in for example _7 , even in the presence of constant magnetic fields and without any periodicity 
condition on V . In the periodic case the situation is somewhat simpler, as we will see in the next 
section. 

Now let us use (|2.7I) and (|2.6p in (|2.16l) and perform formal computations using the completeness 
of Bloch functions and freely interchanging the order of various integrals and series. We would 
arrive at an expression as follows: 



{ H '^Ji J2 ("1 ' k) . . . TTj„n ("n TV 

U,-.u>i (A,, 



-dz 



(k)--z)...(A,„(k)-z) 
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The first question is why is (|2.18p convergent? The answer is not obvious. Moreover, if we look 
at 7rij(a, k) alone (assume without loss that Xj > Xi), we see that using Cauchy's inequality we 
obtain that it could grow like ||(— iQq + fcQ)ui|| ~ y/Xi, which is not very encouraging. For some 
more physical background on such problems, see |14) . 



3 A technical result 

Theorem 3.1. Let M > 1 be an integer and assume that V € C^^^~^(U.'^) is a fl-periodic potential 
with V > 1. Then for every integer < N < M there exists a constant Cn such that uniformly 
m k e 0.* , a e {1, . . . , d} and s, t > 1 we have the estimate: 

|^,t(a,k)| = |7ri,(a,k)| < Cn^tj^- (3.1) 

Proof. During this proof we drop the k dependence of the various quantities. We denote —ida 
with Pa- The eigenfunctions Uj's are a-priori in H'^{VLp) and fJ-pcriodic functions. The lower 
bound on V insures that h>\. The idea of the proof is based on the following identity: 

T^tsia)^ = {h^ut, (pa + ka)h~^Us)T:ts{a) + {ut, h^[pa, h'^]us), (3.2) 

assuming for the moment that the commutator on the right maps into the domain of . 
We know that Pah^^^^ is bounded, thus we can write: 

\7:ts{a)\<CXl/^. (3.3) 

Hence (|3.ip would follow if we can prove the estimate 

\\h^[Pa,h-^]\\<CN. 



3.1 An induction argument 

Lemma 3.2. Choose any cf) £ C^*^(J7p). Then the operators h'^[pa,h~'^] and h'^(j>h^^ are 
bounded for any 1 < N < M . 

Proof. We give a proof by induction in M. Let us start with N = M = 1. Then we have 

hcjjh-'^ ^4> + {-2i{V<t)) • (p + k) - (A(?i)}/i"i 

and 

h[p^,h-^] - [h,p^]h-^ = i{d^V)h-^ (3.4) 

and both are clearly bounded. Note the important thing that we only lost two derivatives of cj) in 
the first relation. 

Now take M > 2. The induction assumption is the following: for any N < M — 1, and for any 
function 4> G C^*^"^(r2p) we have 

\\h''[pc.,h-^]\\<CN, <C^^^, 1<7V<M-L (3.5) 

Note that the induction hypothesis allows us to conclude that 

= /i^0/i-^-(/), iV<M-l, (3.6) 

is also bounded. 

We now have to prove that we can allow N = M. The main idea is to rewrite both h^^[pa, h^^^] 
and h^^ (f>h^^'^ as linear combinations of similar terms but of lower order. We are compelled to 
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perform the induction step simultaneously, because the reduction step mixes the two types of 
operators. 

We start with an identity: 

= h-^[p^,h-^'+^] + [p^,h-^]h-^'+^ = h-^[p^,h~^^+^]+ih-\dc,V)h'^', (3.7) 

where in the second inequality we used p.4p . Multiplying on the left by h^^ we obtain: 

h^'[Po,,h-^'] = h^'~^[pc., h-''+^] + ih^'-\daV)h-'' , (3.8) 

Here the main point is that daV belongs to C^^^~^{flp), thus h^'^^^ {daV)h^^'^^^ must be bounded 
according to p.Sp . 

Now let us deal with the other type of operator. We write another identity: 

= /i-*^</) + /i-*^+i[<^,/^-i] - [h-\[^,h-^'+% (3.9) 

The first three operators in the last equality have their images in the domain of /i*^ due to the 
induction hypothesis p.6p . For the fourth term we apply the Jacobi identity for the double 
commutator and obtain: 

= h'^[h'^^+^,[h,c^]]h-^. (3.10) 

The commutator [h, 0] can be easily put in the form fpPP + 5j with /^'s proportional with 

dp4>, while g is containing up to second order partial derivatives of 4>. Thus: 

/3=1 

+ h-^' {h^'-^[h-^'+\g]h~^) . (3.11) 

The induction hypothesis will insure that the operator inside the large parenthesis is bounded, 
because our /^'s and g belong in particular to (7^^"^. Thus according to (13. 5p and p. 61) we have 
that h^'^~^ fi3h~'^~^^ and h^^~-^[h~'^^~^^ , g] are bounded operators. This concludes the proof. □ 

Remark. Let us show that p.ip is sharp in the free case {V = 0) and for t = s. (Choose the 
lattice to be Z'^). We have (due to the Feynman-Hellman identity) that 27rtt(a,k) = daXtO^) 
at any point k where At(k) is non-degenerate, and At — (27rm + k)^ for some m e Z'^. Then 
|5aAt(k)| ~ |m| - ^/Tt for large t. 

The theorem is definitely not sharp when V — and t s, because in that case all off-diagonal 
matrix elements are identically zero. □ 



4 Applications 

In this section we solve the problems outlined in the introduction, and start by considering traces 
per unit volume. In order to simplify the arguments, we choose V € C°°(f2p). 

Assume that A(k) is a family of bounded, integral operators on L'^{fl), with a jointly continuous 
integral kernel ^(x, x';k), kernel which can be extended to a jointly continuous and periodic 
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function on R^''. The operator B := U* J®. A{k)dk U defined on L'^(R'^) has an integral kernel 
S(x,x') given by (see e.g. [2J): 

S(x, x') = ^ J^^ e*'' (--')^(x, x'; k)dk. (4.1) 

We are interested in a i? operator as given in (I2.17p . Its corresponding A{-) fiber family is: 

A{k) - ^ dzfpD{z){Pc., + fcaj(/i(k) - . . . (p„„ + K^){h{k) - z)-\ (4.2) 

Now if we can prove that Aik) has a jointly continuous integral kernel, then 

Ia,...a„^j:^ [ /^(x,x;k)dxdk. (4.3) 
(^'^) Jn* Jn 

We can integrate by parts M times with respect to z in (j4.2l) using an exponentially decaying 
antiderivative /m of fpD, and write A(k) as a sum of operators of the following type: 



^mi...m„(k) := y^rfz/M(;^)(Pai+A:oJ(/i(k)-z)-"^..(p„„+A:„J(/i(k)-z)-"", (4.4) 
mi + • • • + m„ = n + M. 



Without loss, we may assume that mi > M/n, and keep in mind that M will be chosen to be 
a large number. The integrand becomes trace class if toi — 1 is larger than some number m^ > 1 
depending on the dimension d. The integral kernel of ^,„j^...m„ (k) can be written in distributional 
sense as: 

^mi...m„(x,x';k) (4.5) 

:= X ('^^ ')(A.(k)-.)- ■ ■ ■ (A— (kF^"-- ' 

If we can prove that the above multiple series and the integral are absolutely convergent uniformly 
in x, x' and k, then the integral kernel would be a continuous function as a uniform limit of 
continuous functions. 

Due to standard results on elliptic regularity and Sobolev embedding we can prove the existence 
of m^ > 1 such that 

|u,(x,k)| < C||Mk)'^'^^,(-,k)|| = CA,(k)*^ (4.6) 

uniformly in j and k. Thus we can bound (|4.5p with: 

|^„,...„„(x,x';k)| (4.7) 

<C I d7\U,{7\\ V A rkV^'^'A HcV"''" I^JU2(ai:k)| |7rj-„j„+i(a„,k)| 
<C ^^dz\HAz)\ V,,(k) |,^.^(k)-z|----|A,„^^(k)_,|™„' 

and we will now prove that the above quantity is finite. First, note that l/|Aj(k) — z;| is uniformly 
bounded in j > 1 and z G F. Second, using the triangle inequality and the previous remark we 
have 

iA,(k)-zr-^"^^ + i"i ^ 

again uniformly in j > 1 and z € F. Third, using the exponential decay of \]m\ and the previous 
estimate we have: 

'^''^^"^'^~(kF7F - + 1^1)"'- (4.8) 
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Replacing this estimate in the integral we obtain: 

|^„,...,„„(x,x';k)| (4.9) 
<C A,,(k)-™U,,(k)*'^A,„^,(k)"-|7r,,,,(ai,k)| . . . |7r,„,„^,(a„,k)|. 



As mentioned before, there exists nid large enough such that J2j>i Aj(k)~™'' < oo uniformly in k. 
Now it is time to use p.ip in a certain recursive way. First, we get rid of ji by using the estimate: 

\-nnia^M<C ij^y,^,^^, ■ (4.10) 

Replacing this in (j4.9l) will provide us with a convergent series in ji , the price being the apparition 
of a higher power in Xj^ . 

We can now get rid of all other indices different from j2, starting with jn+i and going down 
to j2- For instance we write: 

l^j„j„+iK:k)| < C f ^ , ^- , (4.11) 

which replaced in (|4.9p allows us to sum over jn+i and obtain the estimate: 

|^^,...™„(x,x';k)| <C ^ A,,(k)-™i+™''+™''+i/2...A,„(k)«"''+i/2|^,_^,Ja„,k)|. 



(4.12) 



Then we can get rid of j„ by writing: 



|^,,._,,„(a„,k)|<C^^^i^i^^^^, n > 2mrf + m, + 1/2, (4.13) 
and then summing. Next comes jn-i- 

X ('k)'~2+i/2 

\nj„_,j„_Aan-iM<C '^^y^ , 7-2 > + n + 1/2, (4.14) 

and so on. Thus we can recursively go down to j2, generating a positive but finite power of Xj^. 
In the end we make use of mi (in fact M) in order to make the series in j2 convergent. 

We have therefore shown that the trace class operator in (|4.4p has a jointly continuous integral 
kernel. Moreover, we now can apply Fubini's theorem and express its trace as: 



dk (4.15) 



^ -.U.(«l,k)...7r,„,,(a„,k) 



Now let us observe that the z integral can be computed through residue calculus, and each 
resulting term will contain a factor of the type (Aj^) with 1 < t < M and 1 < s < n. Each 
such factor decays exponentially. This allows us to prove that each such term converges absolutely, 
by using the exponential decay of /a/ in order to control the growth of the momentum matrix 
elements. As this point we can put together all the z integrals coming from all operators as the 
one in (j4.4l) . Then integrals with respect to z will add up and form the integral 

fFD{z) 



r (A,,(k)-z)...(A,„(k)-z) 
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which win still generate exponentially decaying factors like /]?£)(Aj^), where 1 < t < n and 
1 < s < n. Thus we obtain the final formula announced in (j2.18p . 

Now let us come back to sum rules and derivatives of eigenvalues. While the convergence of 
(I2.12|) is a direct consequence of (13. ip . the absolute convergence of (|2.10p can be proved with the 
same trade-off method as the one we used for the trace problem. 



5 An open problem and an example 

An interesting open problem is to relate the decay of our matrix elements with the regularity of 
V. According to Theorem 13. 11 if V is not infinitely differentiable but has a number of continuous 
derivatives, then one can still prove some decay in energy for the matrix elements. The question 
is whether one can describe the exact decay rate of such a matrix element in case of singular 
potentials. 

Let us investigate a singular "potential" for which one can give an exact asymptotic expression 
for TTij when j ^ oo. Let us consider the one dimensional operator h := —-^ + g5{x)^ g > 0, 
defined through its quadratic form in L^([— i, i]), with periodic boundary conditions. Its domain 
includes the periodic Sobolev space H^{{~^, ^]p) where the respective interval is identified with 
the Id-torus. 

The operator h commutes with the inversion operator. The delta potential does not act on 
the odd subspace, and the odd eigenfunctions are Uj{x) — a/2 sin(2j7rx) , j > 1, corresponding to 
eigenvalues Xj — 47r^j^. On the even subspace, the eigenfunctions must be of the form: 

Uj{x)=Cj{cos{l3jx) + -^sm{Pjx)}, j>l, xG [0,1/2], (5.1) 

where Cj is a normalization constant and Uj{x) Uj{—x) for x € (—1/2,0). Such a function 
cannot belong to iJ^((— 5, \ {0}) unless ii'^ (5-) ^ This provides a quantization condition 
for 0j > 0, which must solve the equation 

/3tan(/3/2)=5, /3 G (27r(j - 1), 27rj), j > 1. (5.2) 

The differences /3j — 27r(j — 1) and Cj — V2 go to zero when j — 00. 

Now let us consider the scalar product ttj := i{ui,u'j) = —i{u'i, uj) and investigate its behavior 
when j grows. Due to symmetry considerations we have: 

■kj — —i1\f2j \ u'lix) siii{2j'Kx)dx = —i m'2^(0+) — i / u'lix) cos{2jTTx)dx 

Jo j"^ Jo 

■ CigV2 . ^/2 fi/^ 
— I — ; h z / (X) fim(2jTrx)ax 



jTT 2j27r2 7o 

■ ^'^"^ + K"(0+) - <(l/2) cos(.,)} + Oil/j 



0{l/f). (5.3) 



..C'i.9\/2 ^ ^^^^ I .3 



Now if we replace in (j2.1ip the operator A with p :— —id/dx^ Um with iii, and u„ with Uj, we 
obtain: 

t{[p,e^'^pe-^'^]u,,u,) = cY,'j^^i^^^+F{t), (5.4) 

where C is a nonzero constant and F{t) is some differentiable function at t = 0. On the right 
hand side we recognize the Riemann function, which in particular is not differentiable sd t — 0. 
Therefore a sum rule like in (j2.12p cannot hold. Nevertheless, the function in (j5.4p is Holder 
continuous at t = of order a g [0, 1/2). 
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